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NTRODUCTION 
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The structure of Input-Output (I-O) modeling has provided the framework 
for organizing data into a set of balanced linear equations (see Leontief, 
1941). The set of linear equations represents the flow of dollars and eneray 
in the U.S. economy. The goal of this modeling exercise is to generate data 
which indicates the energy embodied in a dollar’s worth of goods or services, 
The embodied energy is defined as the Energy Intensity, €, (see Hannan et 
ale; April, 1981), Energy I-O modeling accounts for both the direct and 
indirect energies used to produce goods and service (commodities). Direct 
energy use is de@pinee as the consumption of an energy ebfitiod f0y% such as 
electricity cr coal, at the commodity production site, The indirect use of 
energy occurs, for example, when a sheet of steel is used to produce a car, 
The steel has embodied energy equal to the amount of energy required 
directly and indirectly to produce it. This embodied energy in the steel used 
to manufacture a cay is now accounted for by I-O modeling as an indirect use 


of energy. 


a a ere 


Data are available to model three years, 1943 (sea Herendeen, 1973), 1767 
(see Simpson & Smith, 1975) and 19727 ( see Blazeck et al., 1931). The Bureau 
of Economic Analysis (BEA) is the source of the dollar transaction data for 
the three years. The economic data for 1943 and 1967 were provided at the 
348 sector level and at the 345 sector level for 1972. The three year data 


were formed into a compatible set of 8& sector commodity based constant 
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dollar transaction matrices (see Kennedy et al., July, 1991), 


The data base for each of the three years contains a Make matrix (V), a 
Use matrix (U), and a Final Demand matrix (FD), The V matrix contains 38 
rows which represent 3S industries, Each row contains 35 columns where the 
dollar value of the various commodities produced by the industry are shown, 
The U matrix is also S& by &&, showing the dollar value of the commodities 
purchased by each industry. The label on each one of the rows represents one 
of the Ss defined commodities and each industry forms one column of the 
matrix, The Final Demand matrix shows the dollar value for each of the 38 


commodities purchased by 11 different sectors of the economy. The matrice 


are shown in figure 1. 


cosmodity --> industry --> 
industry commodity 


V matrix U matrix 
88-38 . Bex88 
activity —> 
commodity 
FD matrix 


BBu11 


Figure 1, The Use, Make and Final Demand Matrices Format. 
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Prior to 19772, the dollar transaction data provided by the BEA was 
industry-based, Sectors in the model were defined as industries, ive., by 
establishment rather than by specific commodities produced. This caused 
single establishments to be classified by their secondary products which may 


have been totally unrelated to their primary output. 


The Commodity/Industry assumption, introduced in 1772, explicitly 
acccount for the differences between commodities and industries, The new 
Commodity/Industry model greatly improved the BEAs treatment of secondary 
products by combining analytic and mechanical redistribution techniques (see 


Ritz, i780), 


Revising the 1763 and 1947 data bases was required to make all three 
years compatible. The 1943 and 1947 dollar data were organized into three 
matrices and a vector. The Transfer (TF) matrix showed the dollar value of 
the secondary products produced by industries, The columns af the matrix 
were labeled by industry and the raws by secondary product classifications, 
A Giecet Allocation (DA) matrix indicates the value of industrial output used 
by industries. The columns are labeled by industry and the rows represent 
industrial output (products), The Marginal Direct Transfers (MDT) matrix 
showed the use of marginal or service industries’ output, in rows, by 
industries, in columns. The Value Added (VA) vector is the portion of 


industrial output not accounted for by industry's use of products alone, 


Digitized by the Internet Archive 
in 2023 with funding from 
University of Illinois Urbana-Champaign Alternates 


https://archive.org/details/predictingdirectOOkenn 


The analytical process used to transform the 19763 and 1947 matrices to the 


new form is shown below, 


DA, , + MDT; = Us, (1) 
357 
> Ug cs TSO, (2) 
J=l 
357 

writ VA -— TSO. = TFO. (3) 
vod oS) J J J 


The Total Primary Output (TPO) vector is diagonalized to form a matrix TPO, 


then 
Leone + PERRIS Oe 4 a (4) 
iy 


Several other mechanical date redistribution aperations were performed to 
allow the earlier data to fit the new matrix definitions for U, V, and FD (see 


Kennedy et aly 1951), 


The u and FD matrices are modified ta include the Energy data in BTUs 
(British Thermal Units), The first five rows of U and FD are replaced with 
the Direct Energy Transactions (DET) matrix. The DET indicates the BTU 
Value of each of five energy types used by the industries in the U matrix and 
Final Demand. The data that makes up the DET was found in various 


industrial and gavernmental sources (see Blazeck et aly, Sept. 170), 


| a 4 un ' a 
> +“ ; >i 
i" 7 - 


ae $s > 


© ory : 
en} of esuiviem VdPl bra, Ese) ant milena | 9 babe ee ax 


(r) , Y * 8 hb ot 
s ? eT 3 
eS) ov a ocrF ~ AY ® 

. & t. 


(*3 av = aie + ~~ aT: 
t c Ls . r 


a9 otew anniigyego marhedistadiy 6760 (asinecoge, 


‘rn 7" - 
ee) @% bos (9 eT ro? eeudtinhet spam wen on? ti) of «tenet ae 


ABER yachts 


NTE a wleb veren® att sbulnni of baitthom. em eesirisn ot re us oy 
ifiw betatgs, ete OR fee VU te ewor ovit fest? elt. etiad Lemna it 
OTM off saliotint C20 sat .xtvtem (FS) ehotisaenasT yorend fo 
os xivtem U et) m caleteoiak a? yd beau esqyt yovere svi? — 
euciiay nt baved eke TRO att qu aedem tend steak oat 


Y 


7 i 
ORCL dpe ode fe wresal end) antrwed stat a bya ns 


= 
— 
= he 
ae 
— 
—_ 


The following set of equations shows the mathematical development of the 


I-O model, 


where q; is the row sum vector which is diagonalized to form a matrix q, the 


total commodity output, The matrix is inverted and used to form the Market 


Shares matrix, D; as shown, 


pt (4) 


88 
>, tevedites | ai (7) 
rank s is a 
w= 1 

where g. is the row sum vector which is diaconalized to form a matrix gq, total 


industry output. The matrix is inverted and used to form the Industry Input 


matrix, B, as shown, 


st (8) 


The resulting B and D matrices are in mixed units, The energy commadity 
rows of B are BTU/dollar, while the rest of the matrix is dollar/dollar or 
The D matrix units are dollar/BTU for the energy commodity 


unitless,. 
columns, while the rest of the matrix 1s dollar/dollar or unitless,. 
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The B and D matrices are used in the formulation of the I-Q model equation, 


€ = cr + -ap)7} (9) 
where e is a matrix with one's on the diagonals of the energy rows and 
is otherwise zero. 


The model formulation equation (7) produces the following matrix (see 


figure 2), 


Energy ETU 
Commodities | ET 


Non-Enersy $ 
Commodities | ET 
Figure 2, Energy Intensity Matrix, 

The top five rows of the above matrix, labeled Energy Commodities, are the 
Energy Intensities, Eé, of each of the five energy types Used in each industry, 
The Primary Energy Intensity, eo for each industry is calculated by 
multiplying the Energy Intensity vector for each industry by an A vector (see 
Kennedy et aly, July; 1981), The a, is the fraction of the direct energy for 
the particular fuel type i which is taken directly from the U.S, earth by that 
energy sector. For example, electricity was generated from the primary 
energy types of coal, crude and natural gas, A factor was formulated to 
account for the electricity which was generated from the "ground", such as 
hydro and nuclear power. These quantities were transformed into equivalent 
fossil fuel units at the average efficiency of the existing fossil fuel power 
plants. The resulting foe vector entry represents the portion of electricity 


produced from sources other than primary energy inputs, 
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€é =aA€ (10) 


The final result, eo! is the numerical indicator for the energy cost of 
producing a dollar’s worth of an industry’s primary goods or services, The 
dollar data for 1943 and 1947 has been deflated to 1772 dollar values, so that 
a true interyear comparison of the < can be made, Over time the Energy 
Intensities have changed, as can be noted in the tables in Appendix 1. As the 
price and availability of the 5 energy types have changed due to various 
controllable and uncontrollable events, industry has shifted its use as a 
response. In this way the model thus far could serve as a measure of the 


effectiveness of energy policy and rising eneray prices, 
Mii wey AOI LLY ORrsTHE T-O MODEL 


A desirable feature for any model is that it be able to predict the future, 
In the case of Energy I-O modeling it is desirable that the model be able to 
predict the present, Predicting the present does not appear to be difficult 
until it is understood that the data which is required to predict current 
Eneray Intensities will not be available for 3 to 7 years. From this point of 
view the future is now. The remainder of this document investigates what 


data is required to allow the model to predict the new accurately. 


As a first step the DETs for 1947 and 1972 were substituted for the 1963 
DET in the 19463 model, also the 1972 DET was substituted in the 1747 model, 
The 1943 and 1247 models were then used to calculate the 1947 and 1972Z é, 


These € were then compared to the actual generated for 1967 and 1972, This 
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eperation is referred to as total DET substitution. This is a reasonable 
precedure in view of the fact that complete DET data is available two or 
three years before the completion of the full I/O data for the same year. The 
result of the complete DET substitution is shown on the next three pages in 
histogram form (Figures 2, 3 and 4), Each of the histograms has two plots. 
One plot shows the ability af a complete earlier year model to predict the 
future year, The other plot indicates how well the DET substitution predicts 
the new €, data, For example, the first histogram shows the ability of the 
1743 data to predict the 1947 C, and the ability of the 1943 data with the 
1947 DET substitution to predict the 1947 €, The percent error calculations 
far bath plots on the histograms were performed for each of the #3 sectors as 


follows? 
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HISTOGRAM KEY 
% ERROR DATA FOX 1963 PREDICTION OF 1967 EPSILONS 
Y-AXIS = NUMEER OF EPSILONS SHOWING PERCENT ERROR 
ooo SYMBOL =1963 MODEL PREDICTING 1967 EPSILONS WITH 1947 DET 
Was SYMEOL =1963 MODEL PREDICTING 1967 EPSILONS WITH 1963 DET 


FIGURE 2. 1943 ERKOR DATA HISTOGRAM 
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10 


Z ERROR DATA FOR 1943 PREDICTION OF 1972 EPSILONS 
Y-AXIS = NUNEER OF EPSILONS SHOWING PERCENT ERROR 
ooo SYMEOL =1963 MODEL PREDICTING 1972 EFSILONS WITH 1972 DET 
axe SYHEOL =1963 MODEL PREDICTING 1972 EPSILONS WITH 19743 DET 


FIGURE 32. 1963 ERROR DATA HISTOGRAM 
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Figure 4. 
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abs indicates an absolute value operation 
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In Equations (11) and (12) the superscript number refers to the year of the 
DET data and the subscript number indicates the year of the base dollar data, 
The remaining two histograms, Figures 3 and 4, show data for other year 


combinations as indicated on the plot keys, 


Changing the DET in the I-Q model requires that new H matrix be 
calculated. The B matrix is calculated fram the U matrix with the DET in 
place (see Equation 3) by normalizing the U matrix with the Total Industry 
Output (g). Obviously it is important to normalize the new DET data with a 
new Total Industry Output (g) vector, that represents new total industry 
output for the new energy (DET) data. Obtaining the data to change part of 
the g vector that normalizes the new DET, appears to require little additional 
research, New total industry output data is available from the BEA long 
before the complete transaction matrices are, It was assumed that D matrix, 
which represents an industry’s share of the market, remains constant with 
change in energy use, Therefore, throughtout the remaining calculations the 


D matrix is unchanged from the base year, 
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The total DET substitution experiment proves that substituting complete 


new DET information does improve the current year madel’s ability to predict 


the new Energy Intensity data, 


Table 1 shows the quantitative effect of the 


DET substitution in terms of total absolute value percent error for the 38 


sectors in the model as shown in the following equation, 


The total absolute 


Value percent error is given as a measure of goodness as opposed to mean 


error and 


distribution (see Figures 2, 3, and 4), 


totel Z~ error = 


YEAR TO BE 1947 


PREDICTED 1972 


(% error) 


js. 
i 
— 


YEAR OF DATA USED 

FOR PREDICTION 
1963 1967 

w/o w/ W/O w/ 

1083247920 > a 


13234 752L 10282 6314 


Table 1 


standard deviation because of the skewness of the error 
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Considering that both 1943 and 1947 predicted the following years C, 
equally as well without DET substitution, it is interesting to note that the 
1943 prediction for 19467 with DET substitution produced a better result than 
1967’s prediction of 1972 with DET substitution. This ability could Perce to 
the change in which the 1972 dollar data was collected and the transaction 
matrices formed (see Introduction). Although several mechanical data 
adjustments were performed to make the 1963 and 1967 data compatible to 
1972, the original data from the BEA was gathered under the "industry 


technology" assumption and could still contain inherent differences, 


In view of the positive results found for complete DET substitution, it 
was desirable to determine whether partial DET substitution would provide 
reasonable results. New DET data for partial substitution may be ready for 


use in the model 12 to 24 months before the complete DET can be assembled, 


MOST IMPORTANT DET PARAMETER ANALYSIS 


a a a on a ee 0 i 8 


The Most Important Parameter is essentially a sensitivity analysis of the 
model, If partial parameter substitution is to be sucessful, it seem 
reasonable that the parameters that cause the largest change in the resulting 


should be substituted first, 


The Most Important Parameter Analysis, undertaken here, systematically 
orders the DET parameters in the model according to their effect on 
the €. produced. Previous Most Important Parameters (MIP) studies (see 
Sebald, Nov.,i974, and see Ericson, Sept.,1974) have determined the most 


important parameters for the dollar, as well as the RTU data. However, due 
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to the nature of this research, the importance ranking of the BTU data in the 


DET was done only, 


The method of determining the most important DET involves perturbing 
each DET entry by a fixed percent From the equilibrium state one at a time 
and noting the effect on the € » The percent perturbation was choosen to 
earresond to the average percent change in energy use during the periad 1767 
to 1972, The ehoice of percent perturbation is important because of the 
nonlinear nature of the I/O model, It is possible that another choice of 
percent perce ston could have given different Mast Important Parameter 


ranking, however this phenomena was not investigated, 


In order to determine the ranking of the DET entries, each entry was 
perturbed +30%, The inverse matrix operation shown in Equation ¥ was 
performed to obtain a perturbed set of E,': The perturbed €,’ were then 
compared to the unperturbed set to determine the total absolute value 
percent change caused by the perturbation. The mast important parameter is 


the DET entry that caused the largest total percent change in the €. data. 


1= 


88 | 
MIF % ranking nuaber = 5 (es (€.; = Ei )/e, 100% (14) 


The following table lists the ZO most important parameters, However; a 


more complete list of MIPs is in Appendix Z, 
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When comparing the MIP analysis ranking to the magnitude of the 
parameter entry in the DET, there appears to be no direct correspondence, 
However, a few of the largest magnitude DET entries do appear near the top 


of the MIP ranking, 
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MOST IMPORTANT DET ENTRIES (+30% PERTUREATION) 


COMMODITIES 
PRODUCED 
refined petro 
gas utilities 

electricity 
steel prod, 
Plectricity 
ches prod. 

steel prod, 

ges utilities 
nonferrous prod, 
air trans, 

chem prod, 
crude petro 
electricity 
refined petro 
steel prod. 
electricity 
railroad 
nonferrous prod, 
refined petro 


motor fot. 


Top 20 Most Important Fuel Uses 


FUEL TYPE 

USE FERTUREED 
crude petro 
crude petro 
coz] 

coz] 

Q35 

gas 

O35 

3s 
electricity 
refined 
electricity 
crude petro 
electricity 
refined petro 
refined petro 
refined petro 
refined petro 
gas 

g35 

refined petro 


Table 2 


TOTAL FERCENT 
CHANGE (Equation 14) 
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EFFECTS OF PARTIAL DET SUBSTITUTION 18 


The testing of partial DET substitution was implemented by substituting 
the new most important DET data ane entry at a time and recalculating an 
entire set of Eo This perturbed set Ee, was then compared to the actual & 
and a percent error calculated, The result of substituting the most important 
parameters in order, one at time, for the first 125 important DET entries is 
presented in a graph on the next 3 pages (Figure 5), The graph was expected 
to show that as new DET data was added to the matrix, the resulting Ce 
approximation would steadily became more correct, approaching the error 
level of the full DET substitution, This did not occur, When the third most 
important parameter, the eleetric industry’s use of coal, was substituted, a 
large increase in the tetal percent error occurred (see Figure 5), The reason 
for the large increase is not immediately apparent , especially when we 
eensidered the change in the magnitude of the DET parameter substituted. 
The substitution of this parameter represents only a 15% increase in the 
electric industry’s use of coal. This percent change is average compared to 
ether DET changes. Apparently the single parameter change which the MIP 
“rere identifies is net the camplete reason for a DET entry to be 


eensidered critical to the accuracy of the model. 
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PLOT KEY 


FIGURE 5 - TOTAL % ERROR AS A FUNCTION OF MIP SUBSTITUTION 
Y-AXIS = TOTAL PERCENT ERROR IN EFSILONS 
+++ SYKEOL = NUMEER OF HOST IMPORTANT PARAMETERS SUBSTITUTED 


wxx SYHEOL = NOT USED 
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FIGURE 5- TOTAL % ERROR AS A FUNCTION OF MIP SUBSTITUTION 
Y-AXIS = TOTAL PERCENT ERROR IN EPSILONS 
+++ SYHEOL = NUMBER OF MOST IMPORTANT PARAMETERS SUBSTITUTED 
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I-O MODEL FEEDBACK LOOP ANALYSIS 


Substituting the first 20 important DET entries in various orders indicated 
that particular entries had a greater effect on the Se than the MIP study 
indicated. This effect forced us to determine the importance of the 
relationship between the DET entries. Investigating this phenomenon lead us 
to draw flow diagrams (see Figure 6) and examine critical loops or paths. A 
critical path being one which indicated a large fuel transfer. A critical path 
always contains at least one MIP entry. In eee cases there are multiple 


MIP entries in a path. 


The MIPs are Energy cammodity inputs to various industries, Using the 
industries of the model as nodes of the flow diagram (see Figure 6), it is 
possible to visualize the relationship of the MIP inputs to the less important 
inputs in the model. The diagram shows the five energy producing industries 
and the three industries that are known to consume large quantities of the 
various energy types. It is important to remember that the MIP analysis was 
done for the DET parameters only. Therefore, the effects of feedback paths 
in the remaining part of the Use matrix were not considered and noneneray 


commodities are not shown, 
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Figure 6 I-O0 Model Flow Diesran 
used for viewing Critical Paths 


The electric industry’s use of electricity was Hentiaied to be the 
thirteenth MIP in the model. However, this "self use” feedback path is more 
important to the accuracy of the results, than the MIP ranking indicates. To 
test the parameter’s effect as a loop, it was substituted by itself after the 
ether first 1? MIPs were substituted, The total error in the €, calculation 
was reduced 17.4%. Other "self use" parameter substitutions showed that all 
the "self use" parameters are important to the accuracy of the model, The 
eeal Laceres use of caal had not been identified as a MIP, but proved to be 
impertant when looking at critical paths. Its substitution reduced the error 


in the €o calculation 4%, when substituted after the first 20 MIPs, 


The effects of unbalancing the sum of the energy inputs to energy 
industries can have a greater effect on the results than a cursory look would 
indicate, Consider, for example, the electric industry’s use of coal. If this 
is the only input that changes, as occurred with substituting MIPs one ata 


time, an increase in the use of coal would increase the energy content of the 
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electricity produced, Now that a unit of electricity contains more BTUs, the 
BTU/BTU energy intensity of the electricty output is increased. When the 
RTU of electricity is used in another energy industry, the BTU/BTU energy 
intensity af that energy type is increased, When industry uses any energy 
cammodity input, the energy intensity of the commodity produced increases, 
Clearly, unbalancing the energy input to an energy industry can produce a 
large error due to interdependence of the models components (commodity 


producers), 


The conclusion of the loop and input balancing analysis is that it is 
important to substitute the complete energy industries’ use of energy portion 
of the DET at one time. If that cannot be done; as a minimum, each eneray 
industry’s use of all energy sheuld be substituted. A key to determine 
whether a partial DET substitution has worsened an en estimation is ta laok 
at the resulting change in the energy intensity of the energy commodities, 


These should remain fairly ceanstant as additional DET substitution is done, 


Substituting energy industries’ use of energy first does not preclude the 
fact that the other MIPs in the first 125 group are required to abtain 
reasonable results, Substituting the energy industries’ use of energy merely 


eliminates and explains the initial increase shown in Figure 6, 


amas) atavangnos siebem af% Ia eng sagebsaint af sud Ve159 AGF 


ITE\UTS at? eimabit yrana vaifong ai Leay Rs perineal jaf 


~h te t Tamaya % : iInniots | vi & of jnetzi(es yle ig? nian a4 hivuae 


oh ot ¥ P UPL IEA 4 fugiie {o.sne its 30 pee 2 — 


Sq 


8 
ip — = 
7 


s 
eyes 
Sb 4 ". 


4 y 


‘ i? a 
exon enisined Viiv ale to atin & Vand won Faoubovg nn 


as 


ebeesatant ef idgiua \ seid art? to Htenatat yorens oT Ve 
= + ‘Ai - = Hy 7) 


. 
, > 


aseu yy subd net . beeaewank 2f egy? wovena jadt Yo od 


brane biboinmos eff to volanatiat yorsed af? «fqn yee ne 
L , ee \ % ; 


+= 
aks yrseubet ysvare ae Of Tegal -yersne ef? entcatedsu § (ae 


mo 


- tele 
ec 


ig 


s ; rs } 4 yi bons oop! s(t? 3a oc teulene? alt 
.. 
; -siscgwbat yotees sielomes ad? azote a? tent 
Wo eians ee ot neh ad Tonnes f4dd V2 vem? one Te 7am 


saend? 3 goasvaw ent dabtutijedua Tag tattrag ¢ ain 
‘ > 
¢ 


vr 
muse att 20 ytienetal yprvene eff Mt egneie Qrisuer? 


5 gay ‘setvteube? yutena parrety 


eviup9yt sxe 5% VSL tec) ed? of @8DM wellge Bie jen a 


sane : sent 1 Gh? prt ited saeue -shhvegt sites 
oe inc ovine ney vr lations oct enipiqed clini ete nua 
e y 7 
‘ <a 
4 A 7 ry 


25 


A MATHEMATICAL EXAMINATION OF MIP ANALYSIS 


MIP analysis performs the following mathematical operation 


j= 1. j=i Pe pe 
j=) 1 =: j 0 


cansidering only the effect of a single parameter change on the initial medel 


ta determine the mast important parameter ranking, 


There are eather metheds ef MIP analysis that could be considered, Far 
example, 
440 88 


ae 


J=440-k. i=l 4 Mgt x, 


where x, presents the new initial state caused by the substitution of the 
previous step’s most important parameter. This precess could indeed cause 
parameters that had previous not been identified as MIP to appear earlier in 
the analysis. Inherent in this process is the recognition that MIP are 
dependent on the current state of the model. Avoiding the large increase in 
the total percent error initially shown in Figure 6, could be a result of using 
this MIP analysis algorithm, However, it is not clear that any additional 
benefit can be derived. There is a lewer limit to the DET substitution 
percent error, which is the full DET substitution. If indeed this MIP 
approach could reduce the number of new parameters required for reasonable 


results, the cost for implementing the algorithm is approximately 200 times 
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more expensive to run than the simpler form, 


CONCLUSIONS 


DET substitution is a valid way to get reasonable early estimates of 
Energy Intensities. The histograms in Figures2, 3 and 4 show that a few e 
do not predict well. There appears to be a correlation between these 
intensities and their rate of change over time, Intensities that have 
exhibited a large change from year to year, such as Iron Mining, Chemical 
products, and Railroads, are harder to accurately predict, Inaccuracies also 
occur in industries which are composed of a large number of small firms, such 
as Agriculture, Forest and Fish Services. The diffuse mature of these 
industries reduces the accuracy of the data collection process, allowing 
larger than normal fluctuations in the industries Use data (see Kennedy, et 


arn ieol): 


The ability of the model to provide reasonable results with partial DET 
substitution is indicative of the damped, closed loop nature of the economy, 
LeOry a a eadae in the production af ane commodity doesn’t cause the Grogs 
National Product to become unstable and tend toward infinite. It is 
imperative that the modeler be alert to large data changes in paths that 
contain MIP. While the MIP analysis is not the complete solution to the 
problem of what parameters in the DET must be updated to produce new €,: it 
is an extremely useful tool. The ability to reduce the collective percent error 
in the Energy Intensities produced from 1027% to 703%, by substitution of the 
first 125 MIPs, speaks well for MIP analysis, It is also important that the 


modeler substitute complete energy industries’ use of energy; as described 
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earlier, 


Further research in this could include looking at different methods of 
generating the MIPs. The algorithm in mentioned in the mathematical 
discussion section, which performs parameter substitution during MIP 
analysis, could serve as a ranking index for a MIP analysis, Using the 
ranking ef the maximum change in one €,, eould also be considered. Another 
area of research would be to investigate the change in the Energy Commodity 
Energy Intensities, as an indicator of the amount af unbalance a partial DET 


substitution has caused, 
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APPENDIX 1, ENERGY INTENSITY FOR 1963, 1947 AND 1972, 


SECTOR 


COAL MINING 
CRUDE FETRO 
REFIKED PETRO 
ELECTRIC UTIL 
BAS UTILITIES 
LIVESTOCK 

AG PROD, 

FOREST PROD. 
AG,FOR;FISH SERV, 
TRON MINING 
NONFER MINING 
STONE CLAY riiNING 
CHEM MINERAL MIN 
NEW CONSTRUCTION 
KAINT, CONSTRUCT 
ORDNANCE 

FOOD 

TABACCO 

FASRIC & MILLS 
TEXTILE GO0DS 
APPAREL 

FAG TEXTILE FROD 
WOOD FRODUCTS 
ROOD CONTAINERS 
H’HOLD FURNITURE 
FURN, FIXTURES 
PAFER PRODUCTS 
PAPEREGARD CONT 
PRINTING 

CHEM PRODUCTS 
PLASTICS 

DRUGS 


RUBEER PRODUCTS 
LEATHER PRODUCTS 
FOOTKEAR 

GLASS FRODUCTS 
CLAY STCNE FROD. 
PRIM IKON MAN. 
FRIM NONFER MAN, 
METAL CONTAINERS 
HEATING, PLUMBING 


1963 ENERGY INTENSITIES 


COAL = ERUDE 
1,0045 90062 0033 
60028) = -1,8539 56033 
901320 1,0725 = 11,0830 
1.9710 11,1540 +2820 
80057 14,1005 0098 
8021 30817 21835 
16225 44054 28501 
3071 17139 13582 
8762 27330 17735 
31907 64855 32678 
34143 02609 14520 
24673 38379 27385 
44648 = 202738 41111 
15200 31153 17458 
8920 34184 20145 
15081 1943 B157 
13255 33476 18518 
6288 14893 8419 
32377 62627 26644 
33348 72791 30377 
16005 33025 15049 
24030 AN175 20503 
9098 23403 12934 
17768 20646 13205 
17514 27606 14121 
24402 30458 12346 
u/ F972 83833 304162 
31043 31420 21992 
14895 27916 12057 
94936 = 271867 64099 
73845 197706 71807 
23346 72777 36502 
30370 94793 44205 
5947) 396108 = 372045 
27180) = 317839 =. 284713 
31923 62538 23897 
197101 40323 28539 
12371 24673 41972 
20156 82449 14698 
37862 81672 23132 
116286 70799 27099 
40034 72258 22362 
52914 42615 18818 
42728 41244 16298 


9019 
0011 
60049 
1.1162 
0018 
271 
3207 
B42 
2133 
9724 
12720 
7206 
14974 
2829 
2004 
3684 
3004 
1418 
8887 
B472 
4482 
6545 
2820 
302? 
3879 
3710 
8142 
o447 
3158 
2840 
16433 
3724 
7494 
3433 
9607 
7428 
3444 
2718 
6114 
6778 
7852 
16841 
3173 
6078 
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REFINED ELECTRIC = GAS 


0028 
0043 
0590 
+7022 
1.1109 
8420 
12269 
3632 
8516 
31918 
33425 
27025 
163397 
13785 
944 
16961 
14410 
03/9 
27622 
31660 
14259 
20674 
10027 
12304 
13617 
17196 
45959 
27033 
13775 
124185 
79834 
27289 
35232 
47183 
46986 
28379 
17272 
10721 
63723 
08454 
43495 
49500 
23792 
24477 


PRIMARY 


1.0141 
1.0578 
1.1497 
3.8181 
1.1684 
41849 
08282 
22110 
33706 
125727 
96613 
90281 
267302 
49632 
46116 
37302 
30199 
22728 
102740 
114261 
03137 
77138 
30028 
46778 
30803 
voBI0 
150811 
88231 
45785 
371448 
291203 
102810 
133889 
437607 
366097 
101449 
63404 
37913 
110412 
129770 
193326 
125473 
100707 
87801 
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JOCODE 


4100 
4200 
4300 
4400 
4500 
4400 
4700 
4800 
4700 
2000 
5100 
u200 
3300 
2400 
Du00 
2600 
0/00 
0300 
u700 
6000 
6100 
6200 
6300 
6400 
6501 
6502 
6503 
6504 
6505 
6506 
6507 
6600 
6700 
6893 
6700 
7900 
7100 
7200 
7300 
7300 
7600 
7700 
7800 
7900 


SECTOR 


SCREW MACH FROD. 
FAB METAL PROD, 
ENGINES, TURBINES 
FARM MACHINERY 
CONST, MINING ED 
KAT HANDLING EQ 
HETALWORKING EQ 
SPEC IND MACH 
GEN IND KACH 
KACH SHOP PROD 
OFC, COMFUTER 
SERVICE IND HACH 
ELEC IND AFPARAT 
H’HOLD APPLIANCES 
ELEC LICHT £9 
R-TV COMMUN EQ 
ELECTRONIC COMP 
ELECTRICAL EG 
KOTOR VEH & ED 
AIRCRAFT & PARTS 
TRANSPORT EQUIP 
PROF SCIENCE SUP 
OFTICAL SUPPLIES 
MISC HANUFACTURE 
RATLROADS 

LOCAL TRANSFORT 
MOTOR FET TRANSP 
HATER TRANSPORT 
AIR TRANSPORT 
PIPE LINE TRANS 
TRANSE SERVICES 
COMMUNICATIONS 
R-TV ERDADCAST 
HATER, SANIT SER 
KHOLE,RETAIL TR 
FINANCE SER 

REAL ESTATE 
HOTEL,PERS SER 
BUSINESS SERVICE 
AUTO REFAIR 
ANUSEHENTS 
HED, EDUC SER 
FED GINT ENTERF 
ST,LOCAL GOVT ENT 


APPENDIX 1, CONTINUED 
1963 ENERGY INTENSITIES 


COAL 


34765 
35349 
28711 
34175 
27117 
27073 
22021 
27810 
27087 
19685 
13151 
26481 
24666 
33374 
23826 
13341 
24193 
20527 
31404 
14818 
32199 
15796 
27386 
195 
9937 
£484 
2711 
10044 
4978 
12176 
3472 


11353 


34627 
41134 
31649 
32746 
29236 
27264 
24704 
24961 
29445 
22755 
19709 
34333 
31328 
41048 
35024 
22016 
47949 
31216 
31883 
20643 
37940 
23548 
41412 
35059 
60808 
4190 
31983 

139806 

128032 

150599 
22192 
13228 
10409 
51412 
17884 
12756 

8744 
20830 
16247 
15401 
13269 
15534 
12139 
20412 


14025 
15502 
15905 
13345 
19171 
11426 
10676 
10704 
17188 
9507 
9200 
13957 
11815 
15348 
13270 
10178 
17379 
10895 
13105 
9304 
14510 
9775 
16112 
16088 
55296 
47914 
28476 
134189 
124487 
37975 
17789 
8329 
5196 
31387 
11233 
6178 
5513 
11859 
8044 
849 
4698 
8003 
8254 
12945 


ELECTRIC GAS 
4791-19948 
5590 23398 
4310 15843 
4021 18982 
4168 17032 
3870-15126 
3494 13781 
2620 13860 
4310 16964 
3657 13098 
3014 10036 
5222 19175 
5503 18555 
4162-23924 
4939 20510 
3600 11159 
6696 23294 
5499 18225 
4607 18029 
3868 11079 
451 17979 
3685 12784 
4686-191 
4367-14984 
1467. B17 
287407747 
10734799 
2000 © 42056 
i549 «11074 
5795 64917 
1622 5462 
1685 5310 
1657-5242 
4018 20738 
1925 7118 
1974 6781 
797-3251 
1954 8875 
1984 7792 
1669 6855 
1685 6089 
2557 7697 
906 4301 
4921 14905 
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171716 
27827 
18410 
15527 
71486 
24874 
16875 
12551 
28772 
24381 
23795 
19204 
23247 
15752 
44042 
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JOCODE 


700 


SECTOR 


COAL MINING 
CRUDE FETRO 
REFINED PETRO 
ELECTRIC UTIL 
GAS UTILITIES 
LIVESTOCK 

AG PROD, 
FOREST PROD, 


AG,FOR,FISH SERV, 


TRON MINING 
NONFER MINING 
STONE CLAY MINING 
CHEM MINERAL MIN 
NER CONSTRUCTION 
MAINT, CONSTRUCT 
ORDNANCE 


TEXTILE GOODS 
APPAREL 

FAB TEXTILE PROD 
ROOD PRODUCTS 
HOOD CONTAINERS 
H’HOLD FURNITURE 
FURN, FIXTURES 
PAPER PRODUCTS 
PAPERBOARD CONT 
PRINTING 

CHEM PRODUCTS 
PLASTICS 

DRUGS 

PAINTS 

PAVING 

ASPHALT 

RUBEER PRODUCTS 
LEATHER FRODUCTS 
FOOTHEAR 

GLASS PRODUCTS 
CLAY STONE FROO. 
PRIM IRON MAN. 
PRIM NONFER MAN. 
METAL CONTAINERS 
HEATING, PLUAGING 


APPENDIX 1. CONTINUED 
1767 ENERGY INTENSITIES 


COAL 


1,0048 
0034 
9128 

1.925 
0060 

7247 
7344 
= i) ba] 
i 
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CRUDE 


9055 
1.0515 
1,0417 
1.1191 
1,0730 

37157 

00677 

26877 

20914 

74725 

wasZ8 

w7 376 
188585 

34072 

27405 

26875 

Su616 

18793 

60083 

68724 

32438 

44242 

30796 

25858 

27005 
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84521 

wAS0? 

27147 
2679713 
177416 

u2774 

87928 
377088 
306590 

T8577 

37316 
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82772 
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0925 
10024 
1,0700 
13943 
69049 
26146 
28521 
22596 
11481 
18728 
18727 
24671 
21371 
20384 
18149 
10114 
17895 
7452 
29563 
23715 
12245 
15370 
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12232 
11177 
9474 
34298 
21073 
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SauU? 
01648 
14233 
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274035 
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60013 
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1.1096 
0626 
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3074 
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13061 
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2798 
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4887 
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5038 
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SPY 
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0879 
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GAS 
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wu7 30 
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33289 
170983 
14516 
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1.0114 
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1,1515 
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28292 
118213 
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128880 
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4100 
4200 
4300 
4400 
4500 
4600 
4700 
4600 
4706 
9000 
5100 
5200 
3300 
7400 
5500 
5500 
0700 
0800 
0750 
6000 
6100 
6200 
6300 
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6501 
6502 
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6600 
6700 
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6700 
7000 
7100 
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7300 
7900 
7400 
7700 
7890 
7500 


SECTOR 


SCREW MACH PROD, 
FAB METAL PROD, 
ENGINES, TURBINES 
FARH MACHINERY 
CONST, MINING EQ 
MAT HANDLING EQ 
HETALHORKING EQ 
SPEC IND MACH 
GEN IND MACH 
MACH SHOP PROD 
OFC, COMPUTER 
SERVICE IND MACH 
ELEC IND AFFARAT 
H’HOLD APPLIANCES 
ELEC LIGHT EQ 
R-TV COMMUN ER 
ELECTRONIC COMP 
ELECTRICAL EQ 
MOTOR VEH & EG 
AIRCRAFT & FARTS 
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PROF SCIENCE SUP 
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KISC MANUFACTURE 
RAILROADS 

LOCAL TRANSPORT 
MOTOR FET TRANSP 
WATER TRANSPORT 
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FIFE LINE TRANS 
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R-TV ERDADCAST 
WATER, SANIT SER 
WHOLE, RETAIL TR 
FINANCE SER 

REAL ESTATE 
HOTEL, FERS SER 
BUSINESS SERVICE 
AUTO REPAIR 
ANUSEMENTS 

MED, EDUC SER 
FED GOVT ENTERP 
ST,LOCAL GOVT ENT 


APPENDIX 1. CONTINUED 
1967 ENERGY INTENSITIES 


COAL 


34721 
31141 
20796 
27662 
27595 
23202 
19635 
19712 
23590 
16925 
12546 
23719 
20518 
26799 
20251 
10844 
17744 
i7ill 
23336 
13527 
20328 
14409 
13178 
16406 
202 
6282 
2997 
12263 
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22019 


2208 
2980 
v9? 
17347 
wi81 
3097 
2484 
8646 
3078 
8218 
4387 
8215 
3162 
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CRUDE 


32587 
38316 
28259 
29597 
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26988 
21176 
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30738 
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18443 
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230 
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18507 
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7932 
12173 
10484 

9382 

9037 
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7771 
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7063 

6890 

8315 
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9301 
11444 

9269 

6875 
10631 
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7051 
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BF56 
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11192 
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137796 
125436 
29819 

1629 
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6041 
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3081 
17332 

7850 
15959 
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B476 
13546 
14364 
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SECTOR 


COAL MINING 
CRUDE FETRO 
REFINED FETRO 
ELECTRIC UTIL 
GAS UTILITIES 
LIVESTOCK 

AG FROD, 

FOREST PROD. 
AG,FOR,FISH SERV. 
IRON MINING 
NONFER MINING 
STONE CLAY MINING 
CHEM MINERAL MIN 
NEW CONSTRUCTION 
KAINT, CONSTRUCT 
ORDNANCE 

FOOD 

TABACTO 

FABRIC & MILLS 
TEXTILE GG0DS 
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FAB TEXTILE FROD 
WOOD FRODUCTS 
kOOD CONTAINERS 
H’HOLD FURNITURE 
FURN, FIXTURES 
PAPER PRODUCTS 
PAPEREDARD CONT 
PRINTING 

CHEM FRODUCTS 
PLASTICS 

DRUGS 

PAINTS 

PAVING 

ASPHALT 
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FOOTKEAR 

GLASS PRODUCTS 
CLAY STONE PROD. 
PRIM IRON MAN. 
PRIM NONFER MAN, 
METAL CONTALNERS 
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APPENDIX 1. CONTINUED 
1972 ENERGY INTENSITIES 


COAL 
1.0101 60104 60065 
60028 = =§=61,0554 60056 
60109 =: 1.0622 = 1.0791 
4.9750 1.4465 +6288 
69946 =: 1.0643 0084 
7297 40276 20292 
8399 61448 30607 
9730 68746 97015 
7024 37880 22946 
41540 92965 26578 
24912 61233 28302 
20348 76335 4TSA2 
32425 = 240195 wll 
11878 42274 27088 
7222 4o079 30207 
13755 22810 10367 
9784 40200 19739 
4703 21548 12245 
23090 77310 32083 
18937 76348 31644 
12045 42221 183565 
15916 91613 22103 
10488 45279 24467 
12381 40250 19484 
12981 35074 14085 
17698 37541 14795 
31552 103118 485546 
18554 64962 30492 
9662 SATA 15965 
37794 = 250001 69320 
Su2u2 = s- 171354 67472 
11630 48779 18081 
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32478 45606 17332 


0027 
0016 
60043 
1.1397 
0027 
3871 
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v12i 
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8291 
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CRUDE REFINED ELECTRIC GAS 
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32726 
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21270 
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SECTOR 


SCREW MACH FROD. 
FAB METAL FROD, 
ENGINES, TURBINES 
FARM MACHINERY 
CONST, MINING EQ 
KAT HANDLING EQ 
METALWORKING EQ 
SPEC IND MACH 
GEN IND KACH 
BACH SHOP PROD 
OFC, COMPUTER 
SERVICE IND MACH 
ELEC IND AFFPARAT 
H’HOLD AFFLIANCES 
ELEC LIGHT EQ 
R-TV COMMUN EQ 
FLECTRONIC COMP 
ELECTRICAL ER 
MOTOR VEH & ER 
AIRCRAFT & FARTS 
TRANSPORT EQUIF 
PROF SCTEXCE SUP 
OFTICAL SUFFLIES 
RISC MANUFACTURE 
RAILROADS 

LOCAL TRANSPORT 
MOTOR FGT TRANSP 
RATER TRANSPORT 
AIR TRANSPORT 
PIFE LINE TRANS 
TRANSP SERVICES 
COMMUNICATIONS 
R-TV ERDADCAST 
HATER, SANIT SER 
WHOLE, RETAIL TR 
FINANCE SER 

REAL ESTATE 
HOTEL, PERS SER 
BUSINESS SERVICE 
AUTO REPAIR 
AMUSEMENTS 
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APPENDIX 1, CONTINUED 
1972 ENERGY INTENSITIES 


COAL 


33848 
28347 
23119 
23532 
23808 
21418 
17959 
17774 
21732 
18092 
10532 
19114 
17874 
19853 
17226 
9701 
12893 
16481 
21130 
11974 
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11428 
12441 
14147 
9180 
9312 
2175 
18725 
3407 
12754 
1744 
2449 
w26 
16872 
3145 
1732 
7325 
3930 
9251 
4406 
vee? 
2716 
6940 
11632 


158529 
124813 
BY077 
8013 
16877 
14778 
36923 
11834 
6765 
25979 
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15709 
16518 
12953 
12926 
11403 
12161 
10643 
12741 
12347 
11514 
10575 
12678 
13207 
11371 
13293 
9704 
13129 
14133 
13157 
10478 
15190 
11569 
12104 
16973 
07157 
29236 
33127 
153270 
122449 
30718 
4788 
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B412 
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7126 
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4902 
0181 
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3499 


34 
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APPENDIX 2 
MIP RANKING LIST 


FERCENT ERROR 
841.69 
813.54 
340.61 
217.53 
153.40 


135.87 
93.16 
83.15 
735+35 


73.84 


73249 
70.99 
69.91 
9692 


78.81 


v1.34 
49.84 
40.80 
40.68 
37 +72 


37423 
36.35 
34496 
34.450 
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APPENDIX 2 
MIP RANKING LIST 


FERCENT ERROR 
23493 
ZaeL7 
22243 
Lied 
19.12 


18.95 
18.44 
17.81 
146.20 
14.16 
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14439 
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iz.Z1i 
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